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Abstract 



^sO ' We investigate the properties of charge-neutral /3-equilibrium cold quark 

matter within the Nambu-Jona-Lasinio model. The calculations are carried 

> 

' out for different ratios of coupling constants characterizing the vector and 

scalar 4-fermion interaction, £ = Gy/Gs- It is shown that for £ < 0.4 matter 
is self-bound and for £ < 0.65 it has a first order phase transition of the liquid- 
gas type. The Gibbs conditions in the mixed phase are applied for the case 
of two chemical potentials associated with the baryon number and electric 
charge. The characteristics of the quark stars are calculated for £ = 0,0.5 
and 1. It is found that the phase transition leads to a strong density variation 



^ ■ at the surface of these stars. For £ = 1 the properties of quark stars show 



behaviors typical for neutron stars. At £ > 0.4 the stars near to the maximum 
mass have a large admixture of strange quarks in their interiors. 
PACS number: 14.65.-q, 26.60.+C, 97.10.-q 



I. INTRODUCTION 

The direct application of QCD at moderate temperatures and nonzero baryon densities 
is not possible at present. Therefore, more simple effective models respecting some basic 
symmetry properties of QCD are commonly used. The Nambu-Jona-Lasinio (NJL) model 
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Hill which is dealing with constituent quarks and respects chiral symmetry is one of the 
most popular model of this kind. In recent years this model has been widely used for de- 
scribing hadron properties (see reviews phase transitions in dense matter []5|-[7| and 
multiparticle bound states |8HT0||. This model has been also applied for studying the Equa- 
tion of State (EoS) of /^-equilibrated matter appropriate for stars jn| , [T2] ]. Unfortunately, 
the repulsive vector interaction, which is very important in dense baryon-rich environment 
flU, was not taken into account there. The main goal of our paper is to study the influence 
of this interaction on the EoS and the structure of quark stars. We assume that these stars 
are composed of pure quark matter, similar to previous studies using a simplified EoS of the 
MIT bag model |13,|14||. We do not consider hybrid stars where the quark matter is matched 



to the hadron matter at low densities. The reason of this is two-fold. First, all matching 
procedures are quite ambiguous and model dependent. Second, our calculations show that 
the interiors of such stars have so high densities that the quark degrees of freedom seem 
to be more appropriate. At least this should be a reasonable first approximation to hybrid 
stars. One more advantage of using quark degrees of freedom is related to the problem of 
strangeness content of the star matter. In recent years many calculations have been done in 
hadronic models where strangeness appears through either hyperon formation [113] or kaon 



condensation ]nj -18|]. Within the NJL model the strangeness degree of freedom is handled 
in a very simple way through the population of matter by strange quarks at sufficiently high 
baryon densities. 

II. 13 EQUILIBRIUM QUARK MATTER 
A. Description of the model 



Below we use the SU(3)-flavour version of the NJL model proposed in Ref. |L9|], but 
including not only the scalar but also the vector 4-fermion interaction. The color singlet 
part of the Lagrangian in the mean field approximation can be written as [|IIJ] (h = c = 1) 



2 



£ = J2 V>/ - m f - loGypf) ipf-^fY, Psf 

f z f 

+ ^-Y.Pf + ^X[p Sf + c e . (i) 

Here i/jf is the field operator of quarks with flavour f — u,d,s and 

p S f = <i/>fi>f>, (2) 
Pf = <^/7o</7 > (3) 

are their scalar and vector densities. Angular brackets in Eqs. denote quantum- 

statistical averaging. Gs, Gy and K in Eq. ([[]) are respectively the coupling constants of 
scalar, vector and flavour-mixing interactions. 

The last term in Eq. (P is the leptonic part of the Lagrangian. Below we take into 
account only electrons and treat them as ideal gas of massless particles. As calculation 
show, the maximum value of the electron chemical potential p e does not exceed significantly 
the muon mass so that the muon admixture should be small. In this approximation 

C e = i) e ipif; e , (4) 

where tf) e is the electron field. 

The constituent quark masses, rrif, are determined from the coupled set of gap equations 

m f = m of - G s psf + 2K Y[ Psf , (5) 

f'¥=f 

where mo/ is the bare (current) mass of quarks with flavour / . The single-particle energies 



of quarks are equal to Ef{jp) = yrnj + p 2 . 

Within this model chiral condensates < ■0/0/ > are given by the scalar densities of 
quarks occupying both the negative and positive energy states. The divergent contributions 
of negative energy states of the Dirac sea are regularized by introducing a 3-momentum 
cut-off 0(A — p) in momentum space integrals. In the case of spatially homogeneous matter 
at zero temperature one arrives at the following expression 
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PSf 



vvrif f dpp 2 vrrif 
' 2^~ J EAp) ~ 4^ 

VF f 



PFf 



PFf, 



(6) 



Here z/ = 6 is the spin-color degeneracy factor, ppf 
of quarks with flavour / and 

„2 



^-p/) 1 / 3 is the Fermi momentum 



Vl + x 2 



x~ y \ + x 2 + 1 

— In , 

2 v^T^ 2 - 1 



(7) 



The model parameters rriof,Gs,K,A can be fixed by reproducing the observed masses of 



7r, K , and 77' mesons as well as the pion decay constant /„.. As shown in Ref. ||19|| , a reasonable 
fit is achieved with the following values: 



m 0u = m od = 5.5 MeV, m 0s = 140.7 MeV, 



G. q = 20.23 GeV" 



K = 155.9 GeV" 



A = 0.6023 GeV. 



(8) 
(9) 



Several attempts have been made to extract the vector coupling constant by fitting the 
nucleon axial charge qa 0] and masses of vector mesons pij . It was established that the ratio 
of the vector and scalar coupling constants £ = GyjGs should be of the order of unity. But 
still there is no agreement between the values used by different authors. As demonstrated in 
Refs. P,|i~0|1, the EoS of baryon-rich matter, in particular, the existence of bound states and 
phase transitions is highly sensitive to £ . Due to uncertainty in the parameter Gy , below 
we present the results for various values of £ within the interval < £ < 1 . 



B. Equation of state 

We assume that the quark matter is in chemical equilibrium with respect to the strong 
and weak interactions. If neutrinos accompanying weak processes escape freely = 0), 
one obtains the following conditions 

= Us = Pu + He , (10) 
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relating the chemical potentials fa of various particle species i = u,d,s,e. These equations 
are automatically satisfied if fa are represented as linear combinations of the baryon (ps), 
strange (/is) and charge (/xq) chemical potentials: 

fa = Bifi B + Sifls + QiPQ , (11) 

where Bi, Si and Qi are the baryon, strangeness and charge quantum numbers of the particle 
species i . From Eqs. ([U]) it is evident that ps = and pq = — /i e . In the considered case 
of zero temperature the chemical potentials can be written in the explicit form: 



Pf = y m f +Plf + G vPf , (12) 
fa~p Fe = (3vr 2 p e ) 1/3 , (13) 

where p e is the number density of electrons. 

The chemical potentials, constituent quark masses and particle densities pi are found by 
simultaneously solving Eqs. (|10l)-(|T3|) under conditions of fixed baryon and charge 

densities, 

Pb = Y1 B iPi = \(pu + Pd + ps) , (14) 
i 6 

PQ = J2 QiPi = o ( 2 Pu -Pi~ Ps) - Pe ■ (15) 
i 6 

Unless stated otherwise, the condition of local charge neutrality pg = is assumed. 

To obtain the EoS of quark matter, one can calculate the energy density e = T 00 directly 
from the Lagrangian (|1|). At zero temperature one arrives at the expression (see Ref. (Kj 
for details) 



/ 



PFf 



-^Ylps f + ^-e^. (16) 

The constant e vac is introduced in order to set the energy density of the physical vacuum 
(pFi = 0) equal to zero. This constant can be expressed through the vacuum values of 



constituent quark masses. The latter are obtained by solving the gap equations (^) for the 
case ppf — . At zero temperature the pressure can be obtained from the energy density 
by using the thermodynamic identity 

P = 5>p,-e. (17) 

i 

To characterize the flavour composition of quark matter, we introduce the strangeness 
fraction parameter 

r. = jgU /' (18) 
3B p u + Pd + ps 

where B and S are the baryon charge and strangeness of matter. 

C. Possibility of a phase transition 

Some results concerning properties of the /3-equilibrated quark matter at zero tempera- 
ture are presented in Figs. |1] and |2|. Fig. |I| shows pressure as a function of baryon density 
for three values of the parameter £, namely for £ = 0, 0.5 and 1. Below we compare our re- 
sults with the predictions of the Hadron Chiral Model (HCM) for the /3-equilibrated 



hadronic matter. Pressure calculated in that model is also shown in Fig. |]. Fig. ^represents 
equilibrium concentrations of different quark flavours in /3-equilibrium matter. Comparison 
of the NJL model predictions for different £ reveals a strong sensitivity to this parameter . 
At £ = one can see two zero-pressure points at nonzero baryon densities. They correspond 
to a local maximum and a local minimum in the energy per baryon ej pB ■ In the considered 
case the local minimum corresponds to a bound state of quark matter. This means that 
finite droplets of such matter can exist in mechanical equilibrium with vacuum even without 
gravitational force. 

At £ ~ 0.5 the pressure curve P = P{pb) does not cross zero, but still contains parts 
unstable with respect to the baryon and charge density fluctuations. This implies the possi- 
bility of a first order phase transition of the liquid-gas type. In the case of isospin-symmetric 
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matter similar phase transitions have been studied in Ref. |T0|,^3j . As demonstrated in Fig. 3 
parameters of the phase transition are rather sensitive to £ . Dashed lines in this figure show 
the states unstable with respect to the decomposition of matter into two (k = 1,2) coexist- 
ing phases with different baryon ) and charge (p'q) densities. The spatially averaged 
baryon density is defined as 

p B = \ P % ) + (1-A)pg } , (19) 

where A is the volume fraction occupied by the denser phase (0 < A < 1) . 

The Gibbs conditions of two-phase equilibrium, relating pressures P^ = P (p$ , Pq ), 
and chemical potentials of the two phases can be written as 

p(D = p(2) ; (2Q) 

mS? = /i 2) , (2i) 

$ = $ . (22) 

In accordance with the general discussion of Ref. [p4| , these conditions can be fulfilled 
simultaneously only at nonzero pn . As a consequence, only the "global" charge neutrality 
condition 

Apg ) + (1-A)pg ) = (23) 

should be satisfied in the mixed phase region. By using Eqs. ([T9|) — (|23|) one can calculate the 
equilibrium pressure in this region as a function of baryon density ps ■ As shown in Fig. |3| 
this pressure increases nearly linearly with ps ■ Deviations from the Maxwell construction 
(P = const) are more visible at larger £ . However, the mixed phase region becomes more 
narrow in this case and disappears completely at £ ~ 0.65 . At £ = 1 no traces of this phase 
transition are present any more and the resulting EoS is rather stiff. 



III. PROPERTIES OF QUARK STARS 



In this section the above derived EoS is used to construct the models of compact stars 
composed of pure quark matter. As will be shown below, the properties of such stars differ 
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significantly depending on whether the matter is self-bound or not. In the case of self- 
bound electrically-neutral matter one would expect the existence of macroscopic objects of 
any size above a certain critical mass determined by the surface effects. But certainly the 
gravitational interaction becomes more and more important when the mass of these objects 
grows. 

The gravitational field is described in a standard way by using a spacetime-dependent 



metric g^x") obeying the Einstein's equations [p5|| . We consider only spherically symmetric 
configurations and neglect all dynamical effects, like oscillation and rotation. Moreover, we 
restrict our considerations to stars without magnetic field and at zero temperature. It is also 
assumed that matter can be treated as an ideal fluid. Under these assumptions the Einstein's 
equations are reduced to the Tolman-Oppenheimer-Volkoff (TOV) equations [2(J. For a 



given EoS and a fixed central baryon density p c the inside solution of the metric g^ u { r ), the 
energy density and pressure profiles can be determined by solving the TOV equations until 
the radius R where the pressure vanishes. The outside solution is given by the Schwarzschild 
metric and depends only on the total gravitational mass M of the star. 

The low-density (e < edrip = 4.3 • 10 11 g/cm 3 ) outer layer of a star ('crust') contains 
mainly nuclei and electrons. To describe this nuclear crust we use the EoS suggested in 
Ref. P. 



A. Mass density curves 

Fig. |] represents the gravitational mass of stars as a function of central baryon density p c . 
One can see that the maximum masses predicted by the NJL model depend strongly on the 
relative strength of the scalar and vector interactions £. At £ = the maximum mass has 
a quite low value M max = 1.23 M & with a high central baryon density p™ ax = 9.2 p , which 



is typical for quark matter stars For higher values of £ the EoS becomes stiffer (see 

Fig. [1]) which results in increasing maximum mass. On the other hand, the central density 
p™ ax of the maximum-mass stars decreases with £. For £ = 1 the maximum-mass star has 



S 



M max = 1.60 M and p™ x = 7.8 p . 

In Fig. [| our results are compared with the predictions for neutron stars obtained 
within the HCM []. The upper curve corresponds to the HCM calculation where hyperons 
are neglected whereas the lower one shows the results with the inclusion of hyperons. One 
can see that at higher £ the NJL results for quark stars become close to the predictions of 
the HCM (with inclusion of hyperons) for neutron stars. 

B. Mass radius relations 

The calculated mass-radius relations are presented in Fig. Again the predictions of 
the NJL model are compared with the results of the HCM. One can see that at £ = 1 the 
value of the minimal radius -R m i n = R (M max ) ~ 11.2 km predicted by the NJL model is close 
to the HCM prediction. Here, again we conclude that at large £ the properties of quark and 
neutron stars are similar. 

At decreasing £ the minimal radius becomes closer to the predictions for quark stars 
made in Refs. [|13|,[^J. For self-bound matter (£ < 0.4) the mass-radius relation changes in 
a qualitative way. In this case the corresponding curves start from the origin and M oc R 3 
at small R (without crust). It is not clear to us whether it is necessary to include the crust 
for these self-bound stars or not. But if we do this following the standard prescriptions 
fj~3| we see that the low-mass stars acquire an extended mantel of crust, which may reach 
hundreds of kilometers. On the other hand for quark stars with high masses (near to the 
maximum mass) the inclusion of the crust leads to a relatively small increase of the radii, 
of about several hundred meters. 



1 The slight differences between the results presented here and in [21] are due to improved numerics. 
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C. Density profiles 



Fig. H represents the baryon density profiles psir) predicted by the NJL model as well as 
by the HCM (with hyperons). In all cases the results are shown for the respective maximum- 
mass stars. The profile of a self-bound quark star at £ = has a step-like behavior at the 
star's surface at radius r ~ 8 km. Here the density jumps from ps — 2.6 po to a much lower 
value corresponding to the crust density, or to zero if the crust is neglected. For £ = 0.5 the 
corresponding EoS has a phase transition at a finite (nonzero) pressure (see Sect. [11 C| ). One 
can see that in this case the density in the mixed-phase region changes quite steeply but 
continuously. For £ = 1 no phase transition occurs and the baryon density decreases much 
smother resulting in a bigger radius. The density profile of a maximum-mass star predicted 
by the HCM starts at a much lower central baryon density, p c ~ 6.5 po , but extends to a 
larger radius. Comparing these results we conclude again that with increasing £ the quark 
stars become more and more similar to neutron stars. 



D. Strangeness content 

In this section we discuss the strangeness content of the quark matter as predicted by the 
NJL model. In Fig. ^ the strangeness fraction r s (see Eq. (|I8"D) is shown as a function of the 
baryon density. One can see that the threshold density for the appearance of strange quarks 
depends significantly on the parameter £ and diminishes with it. At p_e>6po the effect of 
the vector coupling on r s becomes less important. For all £ the strangeness fraction tends to 
1/3 that simply means that dense quark matter is composed of equal numbers of up, down 
and strange quarks. This is clear since the difference in bare masses of light and strange 
quarks becomes unimportant at high densities and the SU(3)-flavour symmetry is effectively 
restored. Similar to neutron stars, stable quark stars should have smaller central densities 
than the one corresponding to the maximum-mass star. This stability condition gives a 
maximum value for the strangeness fraction r™ ax in the star's center, indicated by open dots 
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in Fig. 0. These values differ slightly for different £ values. Indeed, r™ ax increases from 0.25 
to 0.29 when £ changes from to 1. In neutron star models dealing with charge-neutral 
hadronic matter the strangeness fraction due to formation of hyperons fT5|f2~T|| is typically 
smaller than in the quark stars. As one can see in Fig. [7], a much smaller value r™ ax ~ 0.12 , 
is predicted by the HCM. 

It is instructive to calculate the total strangeness content of a star. It is obtained by 
integrating the strange number density p s over the star's volume 

Ns — [ \f-~9 Ps u ° dV = 4 7r ( R r 2 — -j= ^ s dr , (24) 

J J o Jl-2m(r)/r 

where u° is the time component of the 4- velocity of matter and m(r) is the gravitational 
mass at the radius r . The ratio of the total strangeness number N s to the total baryon 
number Ng is shown in Fig. ||| for different £ values. Due to the low central density the 
low-mass stars (M<M Q ) have a negligible strangeness content for all £. For larger masses 
N s /Nb is strongly sensitive to £. For the maximum-mass stars the value of N s /Nb changes 
from 0.07 to 0.39 when £ increases from to 1 . This difference originates from the different 
density dependence of the strangeness fraction r s (see Fig. |^) in combination with the 
baryon density profiles of the maximum-mass stars. The stars with higher £ contain more 
strange particles due to a lower threshold in density. The quark stars with large strangeness 
content, say N s /Nb > 0.3 can be named strange quark stars. 



IV. CONCLUSIONS 

In this paper the NJL model is used to construct the EoS of cold /3-equilibrium quark 
matter and to study the structure of compact stellar objects. It was found that the EoS 
depends strongly on the relative strength of the vector and scalar interaction £ = Gy/Gs- 
When repulsive vector interaction is small (£ < 0.4) the matter is self-bound, i.e. there 
is a zero pressure point at a finite baryon density. At higher values of £ (0.4 < £ < 0.65) 
the EoS has no zero-pressure points but still has a first order phase transition at non- 
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vanishing pressure. The properties of the mixed phase were found by using the Gibbs 
conditions for the case of two independent chemical potentials |M] and assuming global 
charge neutrality. We have calculated the mass-density and mass-radius curves, the density 
profiles and strangeness distribution in such stars. It is found that maximum mass grows 
from 1.23 M Q to 1.6 M when £ increases from to 1. In the case of self-bound matter the 
baryon density at the star's surface varies very rapidly from a high value, corresponding to 
the zero-pressure point in the EoS, to nearly zero density. The radii of such stars could 
be quite small even if the crust is added (the minimum radius in this case is about 7 km). 
At £ ~ 1 the quark star global properties (maximum mass ~ 1.6 M Q , radius ~ 11 km) are 
similar to those for neutron stars. However, the quark stars have much higher strangeness 
content as compared with stars composed of hadronic matter with hyperons. These quark 
stars near their maximum masses can be aptly dubbed strange quark stars. 

In the future we are planing to construct a more realistic EoS when hadronic degrees of 
freedom will be included at low densities. 
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FIGURES 




FIG. 1. Pressure as a function of the baryon density (in units of po = 0.15 fm -3 ) calculated 
within the NJL model for different ratios of vector and scalar coupling constants £ . In the case of 
£ = the matter is self-bound that is signaled by the presence of negative pressures. Dashed-dotted 
line shows the predictions of the HCM (with hyperons). 
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FIG. 2. Quark flavour abundances Pf/pB (f=u, d, s) versus the baryon density ps for two 
different ratios of the vector and scalar coupling constants £ = and 1. 
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Pb/Pq 



FIG. 3. Pressure as a function of the baryon density for £ = and 0.5 when a phase transition 
of the liquid-gas type is predicted by the NJL model. The solid lines between the dots correspond 
to the mixed phase obtained by applying the Gibbs conditions. Dashed lines show the results of 
calculations without the mixed phase. 
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FIG. 4. Gravitational mass of a star M in units of the solar mass M & versus the central baryon 
density p c . The three lower curves show the predictions for quark stars within the NJL model for 
different values of parameter £. The dashed-dotted curves correspond to predictions of the HCM 
with and without hyperons. 
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FIG. 5. Mass-radius curves predicted by the NJL model for three values of the parameter £. 
For comparison, the predictions of the HCM with and without hyperons are also shown. The 
results for self-bound stars (solid lines) are presented for two calculations with and without crust. 



19 




FIG. 6. Baryon density profiles for the maximum-mass quark stars predicted by the NJL model 
for three different values of the parameter £ . The dashed-dotted line corresponds to the HCM 
(with hyperons). 
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FIG. 7. Strangeness fraction r s versus the baryon density pb predicted by the NJL model for 
three different values of the parameter £ . The dashed-dotted line shows the results of the HCM 
(with hyperons). The circles correspond to the values of r s at the center of the stars with maximum 
masses. 
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FIG. 8. The ratio of total number of strange quarks N s to the total baryon number Nb as 
function of the gravitational mass M of the quark stars predicted by the NJL model for three 
values of the parameter £ . 
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